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Abstract 

Objective:Jo analyze the association between — 1082A/G polymorphism in interleukin-10 (IL-10) gene and ischemic stroke 
(IS) risk by meta-analysis. 

Methods: We carried out a systematic electronic search in PublVled, BIOSIS Previews, Science Direct, Chinese National 
Knowledge Infrastructure, Chinese Biomedical Database, Weipu database and WANGFANG Database. Pooled odds ratios 
(ORs) with 95% confidence intervals {95%Cls) were calculated to assess the strength of the association. 

Results:? studies were included. There was no significant association between IL-10 — 1082A/G polymorphism and IS risk 
under all genetic models in overall estimates (A vs. G: OR = 1.23,95%CI = 0.85-1. 79;AA vs. GG: 0R= 1. 01, 95%CI = 0.47-2.19; 
AG vs. GG: OR = 0.76, 95%CI = 0.38-1 .55; AA-hAG vs. GG: OR = 0.89,95%CI = 0.46-1 .73; AA vs. AG-hGG: OR = 1 .39, 95%CI = 0.91 - 
2.13). Similarly, no associations were found in subgroup analysis based on ethnicity and source of controls. However, 
removing the study deviating from Hardy-Weinberg equilibrium (HWE) produced statistically significant associations for 
overall estimates under recessive model(AA VS. AG-I-GG OR 1.58, 95% CI 1.04-2.42) and among Asians in all genetic models 
(A VS.G OR 1.64, 95% CI 1.07-2.53; AA vs. GG 0R1.91, 95% CI 1.31-2.80; AG vs. GG 0R1.44, 95% CI 1.09-1.91; AA-fAG vs. GG 
OR 1.54, 95% CI 1.18-2.01;AA VS. AG-fGG OR 1.79, 95% CI 1.07-3.00). Even after Bonferroni correction, the associations were 
observed still significantly in Asians under the two models (AA vs. GG 0R1.91, 95% CI 1.31-2.80, P = 0.0008; AA-fAG vs. GG 
OR 1.54, 95% CI 1.18-2.01, P = 0.001). 

Conc/usion: TW\s meta-analysis indicates that ILIO -1082 A/G polymorphism is associated with IS susceptibility in Asians 
and the —1082 A allele may increase risk of IS in Asians. Considering the sample size is small and between-study 
heterogeneity is remarkable, more studies with subtle design are warranted in future. 
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Introduction 

Ischemic stroke (IS) is a major cause of adult disability and 
death in the world[l], which is a heterogeneous multifactorial 
disease associated with genetic and environmental factors [2]. 
During the past few years, more and more evidence showed that 
inflammatory molecules and the genetic variation of the genes 
which encoded these inflammatory cytokines might take part in 
the pathogenesis of stroke [3] . Inflammatory mechanisms may not 
only play important roles in the manifestation and development of 
IS, but also may be vulnerable to IS in time via accumulation of 
atherosclerotic disease and maintain of atrial fibrillation [4]. 
Several candidate genes of inflammatory cytokines are implicated 
in the pathogenesis of IS, one of which is interleukin — 10 (IL-10). 

IL-10(Gene ID: 3586) is a multifunctional cytokine with anti- 
inflammatory properties, which has been showed involving in the 



inflammatory process of IS [5] . The human IL- 10 gene is located 
on chromosome lq31-32, in which some polymorphisms have 
been found in the promoter region, such as, — 1082A/G 
(rsl800896), -592C/A (rsl800872) and -829C/T (rsl800871) 
[6]. The -1082A/G (also named as -1087A/G in some studies) 
polymorphism could affect IL-10 production [7]. And it is believed 
that the A/G substitution is relevant to low/high amount of IL-10 
secretion, respectively [7]. 

Emerging studies have reported the associations between — 
1082G/A polymorphism in IL-10 gene and IS 
risk[3, 8, 9, 10,1 1,12,13], but the results are inconclusive. Given 
that a single study may be too underpowered to provide reUable 
conclusion owing to relatively small sample size, we performed this 
meta-analysis to estimate the association between ILIO — 1082A/ 
G polymorphism and IS susceptibility more precisely. 
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Figure 1. Flow diagram of the study selection process. 
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Materials and Methods 

Search Strategy 

This meta-analysis conformed to the Preferred Reporting Items 
for Systematic Reviews and Meta-analyses (PRISMA) criteria 
[14]. We carried out systematic literature searches in PubMed, 
BIOSIS Previews, Science Direct, Chinese Biomedical Database 
(http://sinomed.imicams.ac.cn/index.jsp), Weipu database 
(http://cstj.cqvip.com/), Chinese National Knowledge Infrastruc- 



ture (http://dlib3.cnM.net/kns50/),and WANFANG Database 
(http://g. wanfangdata.com.cn/) up to 10th October 2013 to 
identify relevant studies, using the following key words: ("inter- 
leukin-10" OR "interleuMn 10" OR "IL-10" OR "ILIO") and 
(stroke OR "cerebrovascular accident" OR "cerebral ischemia" 
OR "cerebral infarction") and ("polymorphism" OR "mutation" 
OR "genotype" OR "allele" OR "variation" OR "variant"). In 
addition, hand searching of the references in selected literatures 
and the abstracts presented at relevant conferences were 



Table 1. Characteristics of studies included in the meta-analysis. 





study 


Ethnicity 


Study design 


Control source 


Genotype distribution(case/control) 


HWE 


NOS 










AA 


AG 


GG 


iPI 




Marousi2011[8] 


Greece 


Case-control 


Hospital based 


47/53 


71/71 


27/21 


0.72 


6 


Mun5hi2010[3] 


Indian 


Case-control 


Population based 


92/63 


241/218 


147/189 


0.99 


6 


Tuttolomondo2012[9] 


Caucasian 


Case-control 


Hospital based 


58/20 


14/17 


24/11 


0.07 


6 


Jin 201U11] 


Chinese 


Case-control 


Hospital based 


161/78 


27/12 


1/2 


0.09 


7 


Zhang2007[12] 


Chinses 


Case-control 


Hospital based 


202/120 


2/11 


0/0 


0.61 


6 


Lin 2009[13] 


Chinese 


Case-control 


Hospital based 


1 53/83 


28/32 


0/0 


0.08 


6 


Sultana 2011 [10] 


Indian 


Case-control 


Population based 


154/163 


44/47 


40/16 


0.000 


7 
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Study or Subgroup 



Experimental 
Events Total 



Control 
Events Total Weight 



Odds Ratio 
M-H. Random, 95% CI 



Odds Ratio 
M-H. Random. 95% CI 



Jin, 201 1 


349 


378 


168 


184 


12.7% 


Lin 2009 


334 


362 


198 


230 


14.2% 


Marousi,201 1 


165 


290 


177 


290 


17.3% 


Munshi,2010 


425 


960 


344 


940 


19.1% 


Sultana, 2011 


352 


476 


373 


452 


17.5% 


Tuttolomondo,2012 


130 


192 


57 


96 


14.6% 


Zhang, 2007 


406 


408 


251 


262 


4.7% 


Total (95% CI) 




3066 




2454 


100.0% 


Total events 


2161 




1568 






Heterogeneity: Tau^ = 


0.18; Chi^ = 


34.36, df 


= 6 (P 


< 0.00001); 1= = 


Test for overall effect: 


Z = 1.09 (P 


= 0.28) 









1.15 [0.61, 2.17] 
1.93 [1.13, 3.30] 
0.84 [0.61, 1.17] 
1.38 [1.15, 1.65] 
0.60 [0.44, 0.83] 
1 .43 [0.86, 2.38] 
8.90 [1.96, 40.47] 

1.23 [0.85, 1.79] 



83% 



0.01 0.1 1 10 100 

Favours experimental Favours control 



B 

study or Subgroup 



Case 
Events Total 



Control 
Events Total Weight 



Odds Ratio 
IVI-H. Random. 95% CI 



Odds Ratio 
M-H. Random. 95% CI 



Jin, 2011 


161 


162 


78 


80 


7.4% 


4.13 [0.37, 46.22] 


Lin 2009 


153 


153 


83 


83 




Not estimable 


Marousi,201 1 


47 


74 


53 


74 


22.7% 


0.69 [0.35, 1.38] 


Munshl,2010 


92 


239 


63 


252 


26.0% 


1.88 [1.28, 2.76] 


Sultana, 2011 


154 


194 


163 


179 


23.5% 


0.38 [0.20, 0.70] 


Tuttolomondo,20 1 2 


58 


82 


20 


31 


20.4% 


1.33 [0.55, 3.19] 


Zhang, 2007 


202 


202 


120 


120 




Not estimable 



Total (95% CI) 



1106 819 100.0% 

Total events 867 580 

Heterogeneity: Tau= = 0.55; ChF = 21 .84, df = 4 (P = 0.0002); 1^ = 82% 
Test for overall effect: Z = 0.04 (P = 0.97) 



1.01 [0.47, 2.19] 



0.01 0.1 1 10 100 

Decreased risk Increased riski 



Experimental 



Control 



Odds Ratio 



Study or Subgroup 


Events 


Total Events 


Total 


Weight M 


■H, Random, 95% CI 


Jin, 2011 


27 


28 


12 


14 


6.3% 


4.50 [0.37, 54.54] 


Lin 2009 


28 


28 


32 


32 




Not estimable 


Marousl,201 1 


71 


98 


71 


92 


23.7% 


0.78 [0.40, 1.50] 


Munshi,2010 


241 


388 


218 


407 


28.4% 


1.42 [1.07, 1.89] 


Sultana, 2011 


44 


84 


47 


63 


22.9% 


0.37 [0.18, 0.76] 


Tuttolomondo,2012 


14 


38 


17 


28 


18.6% 


0.38 [0.14, 1.03] 


Zhang, 2007 


2 


2 


11 


11 




Not estimable 


Total (95% CI) 




666 




647 


100.0% 


0.76 [0.38, 1.55] 


Total events 


427 




408 








Heterogeneity: Tau^ = 


0.44; Chi= = 


18.56, df 


= 4 (P 


= 0.0010); 1= = 78% 




Test for overall effect: 


Z = 0.75 (P 


= 0.45) 








Fa 



Odds Ratio 
M-H. Random. 95% CI 



0.01 



0.1 



1 10 100 

Favours control 



Study or Subgroup 



Case 
Events Total 



Control 
Events Total Weight 



Odds Ratio 
M-H. Random. 95% CI 



Odds Ratio 
M-H. Random. 95% CI 



Jin, 2011 


188 


189 


90 


92 


6.0% 


4.18 [0.37, 46.68] 


Lin 2009 


181 


181 


115 


115 




Not estimable 


Marousi,201 1 


118 


145 


124 


145 


23.0% 


0.74 [0.40, 1.38] 


Munshi,2010 


333 


480 


281 


470 


27.6% 


1.52 [1.17, 1.99] 


Sultana, 2011 


198 


238 


210 


226 


23.2% 


0.38 [0.20, 0.70] 


Tuttolomondo,2012 


72 


96 


37 


48 


20.1% 


0.89 [0.39, 2.02] 


Zhang, 2007 


204 


204 


131 


131 




Not estimable 



Total (95% CI) 
Total events 



1227 100.0% 



988 



1533 

1294 

Heterogeneity: Tau^ = 0.40; ChP = 20.41 , df = 4 (P = 0.0004); 1= = 80% 
Test for overall effect: Z = 0.34 (P = 0.73) 



0.89 [0.46, 1 .73] 



0.01 0.1 1 10 100 

Decreased risk Increased risk 



Study or Subgroup 



Case 
Events Total 



Control 
Events Total Weight 



Odds Ratio 
M-H. Random. 95% CI 



Odds Ratio 
M-H. Random. 95% CI 



Jin, 2011 


161 


189 


78 


92 


Lin 2009 


153 


181 


83 


115 


IVlarousi,201 1 


47 


145 


53 


145 


Munshi,2010 


92 


480 


63 


470 


Sultana, 201 1 


154 


238 


163 


226 


Tuttolomondo,2012 


58 


96 


20 


48 


Zhang, 2007 


202 


204 


120 


131 



Total (95% CI) 1533 

Total events 867 
Heterogeneity: Tau^ = 0.23; Chi^ = 24.29, df 
Test for overall effect: Z = 1 .52 (P = 0.13) 



13.4% 1.03 [0.51, 2.07] 

15.2% 2.11 [1.19, 3.74] 

16.4% 0.83 [0.51, 1.35] 

18.3% 1.53 [1.08, 2.17] 

17.7% 0.71 [0.48, 1.05] 

13.3% 2.14(1.06,4.32] 

5.7% 9.26 [2.02, 42.48] 

1.39 [0.91, 2.13] 



1227 100.0% 

580 

6 (P = 0.0005); 1= = 75% 



0.01 0.1 
Decreased risk 



1 10 100 

Increased risk 
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Figure 2. Forest plots for association between IL-10 —1082 A/G polymorphism and IS risk in different genetic models. ((A) Allele 
model(A vs. G); (B) Homozygote model (AA vs. GG); (C) Heterozygote model (AG vs. GG); (D) Dominant model (AA+AG vs. GG); (E) Recessive model 
(AA vs. AG+GG)). 

doi:l 0.1 371/journal.pone.0094631 .g002 



performed for other potential related studies. Languages were 
limited to English and Chinese. 

Selection Criteria 

Studies meeting the following criteria were included:(l) the 
study should (■\aluate the relationship between IL-10 gene — 
1082 A/ G polymorphism and IS risk;(2) the study had to be a case- 
control design; (3) genotype distributions in both cases and controls 
were available for calculating an odds ratio (OR) with 95% 
confidence interval (CI); (4)Computed tomographic(CT) or 
magnetic resonance imaging(MRI) were used to assess the 
diagnosis of IS. The following were exclusion criteria: (1) reviews, 
abstracts or animcd studies; (2)studies were not relevant to IL-10 
gene — 1082A/G polymorphism or IS; (3) studies did not report 
genotype frequencies;(4) the design were based on sibhng pairs or 
family. If studies were repeated or overlapped publications, the 
most complete one was included. If studies did not report detailed 
data, we would get in touch with authors to obtain the relevant 
information. 

Data Extraction 

Two investigators {Jin and Li) reviewed and extracted data 
independendy in accordance with the inclusion criteria. The 
results were compared, and if any disagreement appeared, a third 
investigator (Peng) was invited to evaluate such studies, then the 
discrepancy was resolved by discussion. The following information 
were extracted: the name of first author, year of publication, 
country (ethnicity), diagnostic criteria of IS, study design, sample 
size, allele numbers and genotype distributions in cases and 
controls. 

Quality Assessment 

The quality of included studies were assessed by 2 investigators 
(Jin and Li) independently on the basis of Newcastle-Ottawa Scale 
(NOS)[15] which consisted of three aspects: selection, compara- 
bility, and exposure, and each satisfactory answer received one 
point. Studies with a score equal to or higher than five were 
regarded as of high quality. 

Statistical Analysis 

Pooled ORs with corresponding 95% CIs were calculated to 
evaluate the strength of relationship between IL 10 gene —1082 
A/G polymorphism and IS risk under the following five genetic 
models: the allele model (A vs. G), the homozygote model (AA vs. 
GG), the heterozygote model (AG vs. GG), the dominant model 
(AA+AG \-s. GG), and the recessive model {A.\ vs. AG+GG). Z- 
test was used to assess the significance of the pooled OR, in which 
P<0.05 was considered as statistically significant. The Q;test and 
I^-statistics were employed for evaluating the between-study 
heterogeneity, which was considered as significant when 
0.10 or I^>50%[16]. Then, the overall or pooled OR was 
obtained by a random-effect (DerSimonian-Laird method) [17] or 
a fixed-effect model (Mantel- Haenszel method) [18] in the 
presence (Pq<0.10 or I^>50%) or absence (Pq>0.10 or f< 
50%) of heterogeneity, respectively. Furthermore,to explore the 
sources of heterogeneity, we conducted subgroup analysis based 
on ethnicity, source of controls,respectively. Moreover, Bonferroni 
method, controlling for false positive error rate, was utilized to 



adjust for multiple comparisons. As we performed multiple 
comparisons in this meta-analysis for 25 times, the P value which 
was less than 0.05/25 (0.002) indicated statistical significance after 
Bonferroni correction. To validate the reliability of the results, 
sensitivity analysis was performed though omitting one case- 
control study each time, as well as limiting this meta-analysis to 
studies which were conformed to HWE. HWE of genotype 
distribution in the controls of included studies was conducted using 
an onhne program (http://ihg2.helmholtz-muenchen.de/cgi-bin/ 
hw/hwal.pl), and P<0.05 was considered significantly deviating 
from HWE. Publication bias was evaluated by visual inspection of 
symmetry of Begg's funnel plot and assessment of Egger's test[19] 
(P<0.05 was regarded as representative of statistical significance). 
AU statistical analyses were performed using software RevMan 5. 1 
and STATA 11.0. 

Results 

Study Characteristics 

Figure 1 displayed thi; selection process of this study. A total of 
435 literatures were identified after an initial search. Of these 
studies, the first screening excluded 415 citations based on 
inclusion criteria, leaving 20 articles for further selection. Among 
the remaining 20 literatures, 2 explored other diseases instead of 
1S[20,21]; 4 researched other polymorphisms of IL-10 
gene [22,23,24,25]; 2 reported overlapped data [4,8], then the 
one with more complete information was included [8]; 1 was a 
duplicate study[26];l was a review study[27],4 did not have 
sufficient genotype frequencies [28, 29, 30,31]. Finally, a total of 7 
studies were included in our meta-analysis [3,8,9,10, 
1 1,1 2, 13], consisting of 1533 cases and 1227 controls. 

The detailed characteristics of the included studies are listed in 
Table 1. 

AU included studies were of high rjuahty as the NOS score of 
each one was higher than 5 points and the genot)'pe distributions 
in all controls were consistent with HWE except 1 study [10] .All 
the 7 eligible studies were case-control studies, 2 of them were in a 
population-based design[3,10], the remaining were hospital- 
based[8,9,l 1,12,13]. All studies used stroke-free people as controls 
except 1 study [12] which recruited healthy people as controls. AU 
the cases were recruited from hospitalized patients and had a brain 
CT or MRI to assess the diagnosis of IS. 1 study only involved 
first-ever stroke patients [13], 2 studies included first-ever and 
recurrent strokes [3,8] , others did not describe the detaUed 
information [9, 10,1 1,12]. AU studies did not have age limitations 
for cases. Moreover, ethnic groups in these studies were as 
following: 2 were Caucasians [8,9] ,2 were Indians [3, 10], and 3 
were Chinese [1 1,12,13] (n = 3). 

Quantitative Synthesis: Overall study 

As shown in Table 2 and Figure 2,significant between-study 
heterogeneity appeared under all genetic models for overaU 
analysis,thus, a random-effect model was utilized to calculate the 
pooled estimates. OveraU, no significant relationship between IL- 
10 — 1082 A/ G polymorphism and IS risk was found in aU genetic 
models (A vs. G: OR= 1.23,95%CI = 0.85-1.79; AA vs. GG: 
OR= 1.01,95%CI = 0.47-2.19; AG vs. GG: OR = 0.76,95%CI 
= 0.38-1.55; AA+AG vs. GG: OR = 0.89,95%CI = 0.46-1.73; 
AA vs. AG+GG: OR= 1.39,95%CI = 0.91-2.13). 
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study or SubgroLio 
1.2.1 Asian 
Jin. 201 1 
l_in 2009 
Munshi,201 O 
Ziiang, 2007 
Subtotal (95% CI) 
Xotal events 
Heterogeneity; Tstu" = 
Test for overall effect: 



Case 
Events Total 



Control 
Events Total 



Odds Ratio 
Weight M-H. Random. 95% CI 



Odds Ratio 
M-H. Random. 95% CI 



349 
334 
425 
406 



378 
362 
960 
408 
2108 

1514 

0.10; Chi" = 7.43, d< 
Z. = 2.26 (R = 0.02) 



168 
198 
344 
251 

961 
= 3 CP 



184 23. 

230 27. 

940 43. 

262 6. 

161G 100 



1 % 
1% 

o% 

9% 

o% 



1 .2.2 Caucasians 

Marousi,2011 165 290 

Tuttolomondo,2012 1 30 192 

Subtotal (95% CI) 482 
Total events 295 
Heterogeneity: Tau^ = 0.09; Chi= = 2.96. df 
Test for overall effect: Z = 0.22 (P = 0.82) 



177 

57 



290 56.8% 
96 43.2% 
386 100.0% 



B 



Study or Subgroup 



Events Total 



234 

= 1 (P = 0.09); I* = 66% 



Contr 
Events 



1.15 [0.61 . 2.17] 
1 .93 [1.13, 3.30] 
1 .38 [1.15, 1 .65] 
8.90 [1 .96, 40.47] 
1.64 [1.07, 2.53] 



0.84 [0.61 , 1.17] 
1 .43 [0.86, 2.38] 
1.06 E0.63, 1.78] 



Odds Ratio 
M-H. FiJced. 9S% 



0.01 O.I 
Decressed risk 



10 100 
Increased risk 



Odds Ratio 
iVi-H. Fixed, 9S% 



161 
153 

202 

608 



162 
1 53 
239 
202 
756 



2.2.1 Asian 
Jin. 201 1 
Lin 2009 
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Figure 3. Forest plots for association between IL-10 —1082 A/G polymorphism and IS risic based on ethnicity for studies in Hardy- 
Weinberg equilibrium. ((A) Allele model (A vs. G); (B) Homozygote model (AA vs. GG); (C) Heterozygote model (AG vs. GG); (D) Dominant model 
(AA+AG vs. GG); (E) Recessive model (AA vs. AG+GG)). 
doi:l 0.1 371/journal.pone.0094631 .g003 



Notably, there was no GG genotype in either case or control 
group of both the Lin's study and Zhang's study[12,13], resulting 
in that the two studies did not contribute to the pooled ORs in 

homozygote model, heterozygote model,and dominant model. In 
other words, in the three models, there were only 5 substantially 
valuable studies for overall analysis. The main results of this meta- 
analysis were shown in Table 2 and Figure 2. 

Subgroup analysis 

To explore the sources of heterogeneity, we performed further 
subgroup analyses by ethnicity and source of controls respectively. 
Similarly, there were no significant associations in the subgroup 
analyses, and significant heterogeneity in most of the comparison 
models still existed. Table 2 showed the detailed results. 

Sensitivity Analysis 

Sensitivity analysis was conducted to evaluate the stability of the 
results. After excluding the study deviating from HWE[10], 
statistically significant associations were observed for Asians under 
all the genetic models (A VS.G: OR 1.(54, 95% CI 1.07-2.53, 
^=0.02; AAvs. GG: OR 1.91, 95% CI 1. 3 1-2.80,P= 0.0008; AG 
vs. GG: OR1.44, 95% CI 1.09-1.91,P= 0.01 AA+AG vs. GG: 
OR1.54, 95% CI 1.18-2.01,P=0.001,AA VS. AG+GG: OR 1.79, 
95% CI 1.07-3.00, P=0.03, whereas the corresponding pooled 
ORs of total available studies which meant overall analysis for all 
ethnicities were not materially altered (data not shown), except in 
recessive model(AA VS. AG+GG OR 1.58, 95% CI 1.04-2.42, 
P— 0.03). However, after Bonferronic correction for multiple 
testing, the associations revealed still significantiy in Asian 
subgroup under homozygote and dominant models (AA vs. GG: 
OR 1.91, 95% CI 1.31-2.80, P= 0.0008; AA+AG vs. GG: 
OR1.54, P= 0.001). The results indicated that the homozygote 
A.A. and A allele carric-rs (.4.4+ AG) had nearly a 91% and 54% 
increased risk of IS respectively, when compared with the 
homozygote GG in Asians. In particular, removing the study 
which was deviated from HWE eliminated the heterogeneity in 
subgroup analysis of Asians in homozygote, heterozygote, and 
dominant models, but not in the allele and recessive models and 
the pooling analysis (The results were shown in Table 3 and 
Figure 3.).Moreover, omitting the other 6 eligible studies one by 
one, the corresponding pooled ORs in overall comparisons and 
subgroup analysis were not significandy changed and the 
significant heterogeneity between studies stiU existed (data not 
shown). 

Publication Bias 

No obvious visual asymmetry was observed in Begg's funnel 
plots, and the results of Egger's test revealed no statistical evidence 
for publication bias among studies {P= 0.632 for allele model; 
P= 0.773 for homozygote model; P= 0.384 for heterozygote 
model; P— 0.503 for dominant model; P- 0.202 for recessive 
model) (Table 4 and Figure 4). 

Discussion 

It is now accepted that genetics and environmental factors 
contribute to IS susceptibility and outcome. Similarly, the 
inflammation reaction is also relevant to IS. IL-10 is a potent 
anti-inflammatory cytokine with multiple functions taking part in 



inflammation reaction as well as the development of IS. Recently, 
the associations between IL-10 gene —1082 A/G polymorphism 
and the risk of IS have been intensively investigated, however, the 
results are inconsistent. Therefore, we designed this meta-analysis 
to draw a more precise conclusion for the association between IL- 
10 -1082 A/G polymorphism and IS risk. 

In this meta-analysis, no association of the IL-10 —1082 A/G 
polymorphism with IS risk was found under all comparisons, and 
in subgroup analysis by ethnicity or source of controls. However, 
for Asians, after excluding the study deviating from HWE [10], the 
data indicated that IL-IO A allele was associated with increased 
risk of IS in Asians. The inconsistent outcome between Asians on 
subgroup analysis with overaU estimates pardy caused by genetic 
diversity in different ethnicities. In addition, the different subtypes 
of ischemic stroke may contribute to the conflicting results,for 
stroke in Asians is more often due to intracranial atherosclerosis 
than in other populations. Furthermore, as IS is a multifactorial 
disease, except genetic factors, environmental factors also take 
important parts in IS etiology. Thus, this discrepancy may also 
attribute to other environmental factors, such as difierent 
geographic distribution economic status, climate, lifestyle, diet,and 
so on. Importantly, there is currently no consensus for whether to 
include studies deviating from HWE. But if the r('sults are different 
between including or excluding studies deviating from HWE, it is 
suggested that the analysis without studies departed from HWE 
may be more valid [32]. 

Significant between-study heterogeneity displayed among aU 
comparison models. Considering that the diversity in design, 
difference of ethnicity, sample sizes, and measurement errors may 
contribute to common sources of heterogeneity [33], we conduct- 
ed the subgroup analysis by ethnicity and control of sources trying 
to clarify the sources of heterogeneity. Unluckily, we did not 
effectively eliminate the heterogeneity, indicating us that all above 
factors should be taken into consideration. In addition, other 
factors such as subtype of IS, gender distribution, past medical 
history, personal history and so on, might also be responsible for 
the heterogeneity. Notably, after removing the study deviating 
from HWE, the heterogeneity were removed for subgroup analysis 
in Asians under homozygote, heterozygote, and dominant models, 
suggesting that the study deviating from HWE [10] was the main 
source of heterogeneity in the three models for Asians. 

As far as we know, this is the first comprehensive meta-analysis 
exploring the association between IL-10 —1082 A/G polymor- 
phism and IS risk up to now, which involved Caucasian, Indian 
and Chinese populations. In addition, more studies were included 
in our study than a recently published meta-analysis concerning — 
1082A/G polymorphism and IS risk only in South Asians[27]. 
Our meta-analysis also has some advantages. Firstly, the search 
and selection studies were conducted strictiy. Secondly, the results 
of NOS indicated that the included studies were credible. Thirdly, 
no evidence of publication bias was found by Begg's funnel plot 
and Egger's test. Fourthly, multiple testing to adjust for multiple 
comparisons was performed which could reduce the type I error 
rate. In addition, we performed sensitive analysis by excluding 
studies deviating from HWE, considering that deviations from 
HWE in healthy populations may be a sign of selection bias or 
population stratification[34]. 

Despite of the advantages mentioned above, the current study 
has some inevitable limitations that should be acknowledged. First, 
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Figure 4. Begg's funnel plots for IL-10 —1082 A/G polymorphism and IS risk. ((A) Allele model (A vs. G); (B) Homozygote model (AA vs. GG); 
(C) Heterozygote model(AG vs. GG); (D) Dominant model (AA+AG vs. GG); (E) Recessive model (AA vs. AG+GG)). 
doi:1 0.1 371 /journal.pone.0094631 .g004 

there was significant heterogeneity among included studies. Even limiting detailed information such as lacking of subtyping of the 
though we used the random-effects model to calculate pool ORs, ischemic strokes in individual study, we failed to perform further 
the precision of outcome would be affected. Second, owing to subgroup analysis to adjust these possible confounders. Third, the 



Table 4. Publication bias tests for association between IL-10-1082A/G polymorphism and IS. 





Comparisons 


Egger test 




Begg test 




lvalue 


95%Ci 


P value 


A vs. G 


0.632 


(-4.045, 6.047) 


0.548 


AA vs. GG 


0.773 


(-9.428, 7.727) 


0.806 


AG vs. GG 


0.384 


(-7.216, 3.719) 


0.806 


AA+AG vs. GG 


0.503 


(-7.828, 4.812) 


0.806 


GG vs. AG+AA 


0.202 


(-2.163, 7.922) 


0.133 



doi:1 0.1 371 /journal.pone.0094631 .t004 
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sample size in individual included studies was respectively small, 
especially in homozygote, heterozygote and dominant models. 
Forth, only English and Chinese language studies were included in 
this meta-analysis which might have led to bias. Despite of 
conducting an exhaustive search for eligible studies, some relevant 
studies so called "grey hteratures" might be still missed. Fifth, 
similar to a case-control study, meta-analysis was a retrospective 
study, which might lead to recall bias. Last but not the least, this 
study only explored one variation in IL-10 gene, which ignored 
such gene-gene and gene-environmental interactions. 

In [:onclusion,this meta-analysis indicates that IL- 10 gene — 
1082 A/G polymorphism is associated with IS susceptibility in 
Asians and the —1082 A allele may increase risk of IS in Asian 
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populations. However, considering the limitations mentioned 
above, more well designed studies with adequately sized popula- 
tions are needed in future. 
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